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ABSTRACT 


Analytic and experimental investigations of the interactions of 
hig j-energy protons with spacecraft materials and light nucleus atoms 
have been conducted. Irradiation of typical materials for scintillation 
detectors using accelerator-produced energetic protons and associated 
theoretical calculations indicate activity within detectors produced by 
this mechanism is significant to detection of low-level cosmic gamma 
radiation. Samples of several materials were flown aboard Skylab to 
obtain measures of the neutron and proton flux environment. Extremely 
weak activity waB produced in the samples. Calculations of interactions 
of intermediate-energy protons with light nuclei were performed and a 
more complete theoretical model programmed for digital computation. 
Preliminary results appear promising for calculating products of such 
interactions, but additional work remains to be done. 
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1. INTRODUCTION 


This report synopsizea th vaults of a series ot investigations, 
both analytic and experimental, of the effects of the radiation environ- 
ment on Earth-orbital and balloon-borne systems and a study of cosmic 
sources of relatively high-energy photon radiation. 

The radiation environment was considered for its important 
inducing radioactivity in such a manner as to affect the capability's 
of radiation detectors designed to study distant low-level gamma-ray 
sources. This radioactivity can be induced in the structural materials 
of the measurement platform, the atmosphere (for balloon-borne measurement 
systems), or in the detector material itself. 

Both the early and late work performed on this contract has been 
dedicated to theoretical analyses of reaction cross-sections. The 
earlier results were documented previously; this document provides a 
summary of the latest work. 

Another portion of the work waB directed toward a combined 
experimental and computational program of investigation of induced radio- 
activity determination. Radioactivity induced by protons in sodium 
iodide scintillation crystals was calculated and measured directly. These 
data are useful in determining trapped radiation and cosmic-ray-induced 
background counting rates in space-borne detectors. Neutron- and proton- 
induced radioactivity at various locations within Skylab was measured. 
Low-level gamma-ray spectroscopy measurements were made on the returned 
samples to determine the Incident neutron and proton fluxes and spectra 
and their variations with mass distribution. 

Teledyne Brown Engineering also contributed to the design of a 
balloon-borne gamma-ray measurement apparatus and the construction and 
checkout of the apparatus and its use for astronomical measurements. 
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2. EARLY THEORETICAL WORK 


A summary of the status of the problem of describing high-energy 
nuclear reactions Including descriptions of specific models to predict 
observed phenomena waB presented (Ref. 1), This work wqs extended in the 
last few months and formed the basis of the experimental work that fol- 
lowed the original analysis. A more complete description of the analysis 
and recent modeling work are presented in Section 4 of this report. 


Z-\ 


?. INVESTIGATION OF INDUCED RADIOACTIVITY 


3.1 STUDIES 01’ INDUCED RADIOACTIVI' ■ IN DETECTORS 

Proton-induced radioactivity of scintillation detectors iB an 
important source of background counting rate whenever these detectors ore 
used in a space environment. In addition to the ever-present primary 
cosmic-rays penetrating spacecraft in Earth orbit, the trapped proton 
radiation adds considerably to the total proton flux. For orbitB of low 
altitude and medium inclination, the largest fraction of proton flux is 
received in the region of the South Atlantic Anomaly. Figure 3-1 shows 
the flux contours for protons with energies greater than 5 MeV at an 
altitude of 600 kilometers (Ref, 2). Figure 3-2 gives the proton energy 
spectrum averaged over 24 hours for an orbit at an altitude of 720 kilome- 
ters and an inclination of 30 degrees (Ref. 3). Figure 3-3 gives the 
primary cosmic-ray spectrum near the Earth. 

Radioactivity induced in spacecraft scintillation detectors was 
first recognized and analyzed by Peterson on his 0S0 1 experiment (Ref. 4). 
The analysis of data from the Russian gamma-ray experiments on the 
Cosmos 135 and Cosmos 163 spacecraft has included corrections for the 
radioactivity induced by passage through the region of the South Atlantic 
Anomaly, although the details of these corrections have not been given 
(Ref. 5). Dyer and Morfill (Ref. 6) have calculated the induced radio- 
activity in CsI(T2.) by the same general method used by Fishman (Ref. 7) 
and have compared these calculations to radioactivity induced by 155-MeV 
protons. This energy is representative of that of nuclear-interacting 
trapped protons. Their measurements are used in a study of the anticipated 
background radiation on the UK-5 celestial gamma-ray experiment. The 
data contained within Reference 7 have been used to determine the effect 
of cosmic-ray-induced radioactivity on diffuse celestial gamma-ray 
measurements (Ref. 8). 

To deduce the amount of radioactivity and its energy and Lime 
decay spectrum induced in scintillation detectors by protons, both an 
analytical approach and an experimental approach were undertaken. The 
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FIGURE 3-1. PROTON FLUX CONTOURS FOR E D > 5 MeV 
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FIGURE 3-2. PROTON ENERGY SPECTRUM FOR h = 720 km, i = 30 degrees 
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FIGURE 3-3. PRIMARY GALACTIC AND SOLAR COSMIC -RAY SPECTRA 
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analytical approach consisted of calculation# of spallation yields and 
a study of the nuclear properties of the spoliation products. The 
experimental approach waB to irradiate a scintillation detector with 
high-energy protons and directly observe the counts resulting from the 
decay of induced radioactivity. Both approaches are limited in their 
estimation of the radioactivity induced in u true apace environment. 
However, within errors, both estimates were in agreement and may prove 
useful in design and data reduction applications for space-borne experi- 
ments using scintillation detectors. 

Calculations of spallation yields were made for protons inter- 
acting with sodium iodide at energies above 100 MeV using the semi- 
empirical formulas derived by Rudstam (Ref. 9). These formulas give 
the cross section for spallation production within a factor of ihree 
over a wide range of atomic mass numbers and incident proton energies. 

I* 1 i n effort to measure directly the radioactivity induced in a 

scititil'r detector, an accelerator experiment was performed at the 

NASA Space Radiation Effects Laboratory, Newport News, Virginia. A 

sodium iodide detector, 17 gm-cm -2 thick, was Irradiated with a low flux 

of 600-MeV protons over a 10-second period. Tim integrated irradiation 

flux, 7 * 10 1 ® protons in a 5-centime Ler-diameter sodium iodide detector, 

127 123 

was measured using the known cross section for the I (p» p An) I 
reaction as a beam monitor (Ref. 10). The number of 1^3 nuc i e i pro- 
duced were derived from the 0.192-MeV line observed in the decay spectrum 
at a known time after the irradiation. Continuum and background radiation 
corrections were made to the data in this derivation; however, the accuracy 

of the data is estimated at only ±30 percent because of the magnitude 

1.23 

of the corrections and the uncertainty of the I cross section. 

Counting rate spectra were taken at various times after the 
irradiation. In general, the spectra show exponential cuntinua with 
several discrete peaks superimposed. Most of the peaks from electron 
capture and internal transition decays of radioactive isotopes of iodine, 
tellurium, and antimony could be identified in tiie spectra. The continuous 
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nature of the spectra arises from the large number of radioactive 
products involved and the fact that most transitions Include a non- 
discrete beta transition. 

The study indicated that pro ton- induced radioactivity will be 
a major source of background radiation for low-level, space-borne 
detectors in the hard X-ray and gantma-ray regions. Approximately 
one gantma-ray will be produced after each nuclear Interaction. The 
decay of the Induced radioactivity will roughly follow a 1/t rate, 
and the energy spectrum is expected to be exponential with an e-folding 
energy of approximately 0.9 KeV, Proton Irradiation data at 600 MeV 
are in agreement with the above results. 

More precise calculations may be made using experimental cross 
section data and detailed decay schemes. The actual proton transport 
within heavy spacecraft will be needed for these calculations, along 
with the proton dose history. 

3.2 STUDIES OP INDUCED RADIOACTIVITY IN SPACECRAFT MATERIALS 

Large, orbiting spacecraft are expected to produce significant 
numbers of neutrons through nuclear reactions between trapped protons 
and spacecraft materials. Although measurements of upward-moving 
atmospheric neutrons (the neutron albedo) have been made, locally 
produced, secondary neutrons in spacecraft have not yet been studied. 
Secondary neutrons have previously been shown to contribute a negligible 
dose to astronauts, but their great penetrating power and ability to 
activate many materials may seriously degrade the sensitivity of ) -ray 
and gamma-ray detectors on orbiting spacecraft. This is especially 
true of scintillation and solid-state detectors aboard the High-Energy 
Astronomy Observatory spacecraft. For these detectors, prompt, charged- 
part icie-produced background can be greatly reduced by active shielding, 
whereas there is no known method of reducing activation background in 
the detectors. 
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The Skylab Neutron Environment Science Demonstration measured 
neutrons and protons by the activation analysis technique (Ref, 11). In 
this technique, samples of known composition were exposed to the unknown 
neutron and proton environment at various locations within Skylab, After 
return of the samples, they were analyzed in a low-level, high-resolution 
gamma-ray spectrometer to deduce the type and amount of radioactive Isotopes 
produced during their irradiation. The neutron and proton environments 
can then be derived by using the appropriate activation cross sections. 

Four activation sample packets were flown and returned on the 
Skylab XV mission. They were deployed by the crew at the locations 
within the Orbital Workshop (OWS) given in Table 3-1. The four activa- 
tion sample packets each contained five samples of high-purity metals, 
with dimensions of 1.91 by 1.91 by 0.32 centimeters. Tour of the metal 
samples were tantalum, nickel, titanium, and hafnium. The fifth sample 
was tantalum with a covering layer of 0.65-millimeter-thick cadmium, 
which excluded thermal neutrons. The cadmium had a minimum absorption 
length of seven mean-free paths for neutrons with energies less than 
0.2 eV. The samples were sewn into beta cloth (fiberglass) to form 
the activation packets. Table 3-2 indicates the properties of the 
samples . 

The activation packets were in the Skylab IV orbit for the mission 
duration (Novel ,jr 14, 1973, to February 5, 1974). During the initial 
5.3 days and the final 2.7 days of the mission, the packets were stowed 
in a Command Module locker. The remainder of the time, 76.0 days or 
91 percent of the orbital time, the packets were deployed at the locations 
given in Table 3-1. 

The radioactive isotopes found in the returned samples, along 
with their nuclear decay properties and production modes, are listed in 
Table 3-3, All activities found were extremely weak; the average count 
rate for a particular isotope was less than one count per hour. 
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TABLE 3-2. ACTIVATION SAMPLES 


SAMPLE 

NOMINAL 

WEIGHT (grams) 

MAIN 

ISOTOPES 

1. Tantalum (bare) 

19 

Ta 181 (100%) 

2. Tantalum (cadmium-wrapped) 

19 

Tal81 (100%) 

3. Nickel 

12 

NlS8 (68%), Nl60 (26%) 

4. Titanium 

5.6 

T1 48 (74%), T1 46 (8%), 
T1 47 (7%) 

5. Hafnium 

17 

HflBO (35%), Hf] 7 ? (14%), 
Hf 178 ( 27 %), Hf^ 77 (19%) 














TABLE 3-3. ACTIVATION ISOTOPES IDENTIFIED 
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The preliminary measurements showed that both neutron and proton 
activation isotopes were produced in the samples from which several 
qualitative observations can be made: 

58 58 

• Co , which has a high-yield cross section from Ni 

via the fast neutron (n, p) reaction, is produced in 
greater quantity in the higher mass locations, 

• Co'’** and V^, which are proton activation products of 
nickel and titanium, respectively , are produced at the 
lower mass locations. Co^? shows no clear enhancement. 

• Tal82, a neutron capture product, is produced more at 
the higher mass locations. The capture cross section 
rises rapidly at thermal neutron energies. The lack 
of enhancement at the water tank location probably 
indicates a low thermal to fast neutron flux ratio. 
Measurements of the cadmium-wrapped tantalum, not 
made under this contract, should Indicate the thermal 
neutron component. 

« Proton activation in the film vault location, as 

evidenced by Co^b and V^, indicates the presence of 
the high-energy tall of trapped protons and/or primary 
galactic cosmic rays. Some of this activation could 
be due to temporary storage of the packet while in the 
Command Module. 

Reference 11 contains more complete information on this 
experiment. 
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4. RECENT THEORETICAL TREATMENT 


4.1 INTRODUCTI ON 

Charged particles, protons, and alpha particles, and in some cases 
light nuclei such as gC* 2 ami are effective tools to study nuclear 

structure by investigating the reaction products with target nuclei in 
various mass ranges. In particular, proton reactions with selected light 
nuclei such as gO 16 , and ^qNc 20 are considerable interest 

because of the abundance of these nuclei in the upper atmosphere. Reactions 
with these nuclei by protons occurring in cosmic radiation in the upper 
atmosphere are of concern for balloon-borne measurement apparatus. Of 
particular Interest are proton reactions with light nuclei for proton 
energies ranging from 10 to 20 MeV, It is of interest to investigate the 
mechanisms of these types of reactions by developing an appropriate 
formalism and to study the gamma-ray production cross section as a function 
of the incident proton energy. The most recent work on this contract 
was again addressed to the problem of these computations. 

4.2 THEORETICAL MODEL 

4.2.1 Types of Reactions 

When a charged particle such as a proton approaches a nucleus with 
energies low compared with the Coulomb barrier exhibited by the target 
nucleus under consideration, the particle interacts only with the 
Coulomb field, and this leads to ordinary elastic scattering. As the 
energy of the incident proton increases, it begins to interact with 
the nuclear field so that the approaching particle can be either elastically 
scattered or absorbed to form a compound nucleus. If a compound nucleus 
is formed, this intermediate nucleus will be in an excited state from 
which de-excitation will take place either by emitting a proton with 
energy equal to that of the incident proton or by emitting a proton of 
energy different from - usually less than - that of the incident proton. 

The former type of reaction is called a compound elastic reaction, and the 
latter is called a compound inelastic reaction. In the compound inelastic 
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reaction, the difference in energy between the incident and emitted proton 
is left in the residual nucleus, leaving it in an excited state; it 
de-excites through gamma-ray emission until it reaches its ground stjte. 
However, if the incident proton energy is sufficiently large, then the 
compound state can be followed by particle emission leading to a residual 
nucleus different from the target or the compound nucleus. The residual 
nucleus is generally left in an excited state, and thi3 emission of a 
particle different from the projectile particle leads to a de-excitation 
mechanism, which could open various channels of particle and gamma-ray 
emission. If the original reaction leads to emission of more than one 
particle, the mechanism is rather complicated and a simple approach to 
this aspect of the problem is not possible. 

In addition to elastic and compound nucleus rt actions, a third 
type of reaction chat can take place is the direct reaction. In the 
case of direct reactions, the Intermediate compound nuclear state is not 
formed; instead, the projectile reacts directly with the target nucleus, 
leading to inelastic (p,p') scattering or pickup or charge exchange 
reactions. A typical pickup reaction is (p,d), and a similar charge 
exchange reaction is (p,n). These different types of reactions arc 
represented schematically in Figure 4-1. 



NUCLEUS 


RESIDUAL NUCLEI 


COMPOUND 

NUCLEUS 


FIGURE 4-1. COMPOUND NUCLEUS ENTRANCE AND EXIT 
CHANNELS 
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4,2.2 Elastic Scatterin 


The theory of elastic scattering is very well worked out and is 
rather well understood. However, for the sake of completeness, a brief 
outline of the theory is described here. The reaction mechanism of a 
spinless neutral purticle iB discussed first so that the discussion of 
charged particle reactions can ho understood better. The Schrodinger 
equation for such a reaction can be written as 

7^ if* + [E - V(r) ] ij> = 0 . (4-1) 

The solution for this must satisfy, for large r , an incident plane wave 
and a scattered spherical wave. 

Assuming that the particle is approaching the target nucleus in 
the Z direction, one can write the asymptotic form of the wave function 
in spherical coordinates as 

«|/ - exp (ikz) + H P— ( li yA f( 0) f (4-2) 

where 

. /2mir 

k " fi 

m is the mass of the projectile, and E is its energy. The problem is to 
calculate f(0) for an appropriate potential. To do this, one attempts 
to solve Equation 4-1 by expanding the wave function into a product of 
radial and angular parts and writes 





U £ (r) 

(cos 9) 


(4-3) 


where i - 0, 1, 2 ... represents the orbital angular momentum quantum 
number and I 3 £ (cos 0) are ordinary Legendre polynomials. 

When as given in Equation 4-3 is substituted into Equation 4-1, 
Equation 4-1 will lead to 
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r 


U* “ 0. 


(4-4) 


d* U, 
dr 2 


+ 


— (E - V) 
*2 


Z(Z + 1) 
r 


The angular part disappears and it i8 only the radial part that appears 
in Equation 4-4. To make matters simpler, one at this stage looks into 
s-wave scattering, which immediately reduces Equation 4-4 to a simple 
form since & is set equal to zero. Thus, 


d z U 0 
dr z 



(E - V) U 0 - 0 


or 


d 2 U 
dr 2 " 


k 2 U 0 


0 . 


(4-5) 


Obviously, Equation 4-5 lias a solution of the form 
Uq •• A sin (kr + 6) 

«= A exp (-16) [exp (216) exp (ikr) - exp (-ikr)]. (4-6) 

However, from Equation 4-3 for s-wave scattering, one obtains, since 
Pq (cos 8) is equal to unity, the relation 

ID = “ . (4-7) 


Therefore, for Z - 0 


<* # in + 


out 


(4-8) 


where gives the incident, wave and *l' out gives the scattered wave. 
However, by the assumption made in Equation 4-2, the. incident wave w^ll 
be given by 


i|/ in = exp (ikz). 


(4-9) 
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For e-wave scattering, this can be determined from the exp.in n 
of Equation 4-9 into a combination of Bessel functions and Legendre poly- 
nominals and by setting A to zero in the expansion obtained. This leads 
to 

«0 

exp (ikz) - y (2A + 1) i ? ' j*(kr) P* (cos 0), (4-10) 

AO 

so that for s-wave scattering Equation 4-10 can be rewritten, by making 
use of the expansion of Bessel functions and also by remembering Pg (cos 4) 
is unity, as 

exp (ikz) U - 0) - (4-11) 

or 

exp (ikz) (A = 0) * 2ikr ^ exp ^kr) - exp (-ikr)] . (4-12) 

From physical considerations, it is obvious that the scattered wave is 
given by 


4- out (A = 0) = if* - ip in (A = 0) . 


(4-13) 


At this stage, one uses the solution of the Schrodinger equation 
as given in Equation 4-6, the relation given in Equation 4-7, and 
Equation 4-12 to obtain 

4-out (A = o) - A (exp (216) exp (ikr) - exp (-ikr)] 


‘Ijgr [exp (ikr) - exp (-ikr)]. (4-14) 


Again, the physical condition requires that the outgoing wave should con- 
tain no exp (-ikr) term, and this condition can be satisfied in Equation 4-4 
by setting the terms containing exp (-ikr) to add up to zero. Such a 
restriction leads to 
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“§ » !f(0)| 2 . (4-19) 

The value of this Is determined by making use of Equation 4-18. Of course 
this will give only the value for the s-wave scattering. However, using 
a similar approach, one can evaluate the cross section for p-wave (1*1) 
or higher-order scattering cross sections, 

4.2.3 Compound-Nucleus formation and Uauser-Fesbach Formalism 

The discussion in the previous section of shape elastic scattering 
now leads to compound-nucleus reaction as the next problem in sequence 
to be confronted. The simplest type of compound-nucleus reaction is the 
one in which no direct reactions exist in the nonelastic channels, and in 
such a case a statistical approach is valid. As in the case of elastic 
scattering, it is convenient to discuss the reaction from one initial 
state to one final state in the absence of spin. The details of the 
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formalism have been worked out by Hauser nml Fesbuch and also extensively 
dlBcusBed In the literature. Therefore, the salient aspects of the 
theory will be briefly discussed. 

Let a denote the Initial state and & denote the final state , and 
let the spins of the entrance channel be denoted by i and 8 and similar 
quantities for the exit channels by l' and s* . Here, s denotes the 
vector sum of the spin of the projectile and target nucleus and s' denotes 
the same for the residual nucleus and the emitted particle. Thus, one 
can write 

s + l » s' + i' (4*20) 

from conservation of angular momentum. 

The compound nucleus can decay elastically, in which case the 
emitted particle is identical to the projectile. In the terminology 
above, this implies that 

u. « f>. 

The scattering matrix, S, in this case will be diagonal. As before, 
if one begins by considering the reaction in the case of spinless particles, 
the cross section for the formation of the compound nucleus can be shown 
to be 



where Is the scattering matrix. The cross section for the formation 
of a compound state leading to an emission back into the entrance channel, 
which implies compound clastic scattering, may be obtained from Equation 
4-21, where S. 1, leading to 

°aa * 11 O - 1 S ua| 2 J 

■ II T . (4-221 

a a 
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In Equation 4-22, 1b the transmission coefficient. Transmission 
coefficients can be calculated by using nn appropriate optical potential. 

In the nonelastic case, one la Interested in o a g, and to evaluate this it 
is necessary to determine S u „j 2 . The standard technique used is to upply 
the physical principle that the formation and decay of the compound 
state are Independent of each other so that one can write 


<M> • <M> 


(4-23) 


which is equivalent to stating that 


*a u a6 - *B “SB • 


If the probability of decaying into channel ft is designated by 
Pg, then Equation 4-21 can be written as 

CN 


P 6 


JI ** T u Pg . 


(4-24) 


Following the theoretical formalism of Hauser and Fesbach, the probability 
of the formation of a compound nucleus of spin J from channel a is 
(from Equations 4-22, 4-23, and 4-24): 


T a p g “ T a P p 


(4-25) 


lliis statement is true for all possible channels u and S. One can con- 
clude, from Equation 4-25, that 




T A 




(4-26) 


where X is a constant. Obviously, the probabilities hold the relationship 


I - I 


p «* i 

ft 
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and 


I P c - * 1 T « - *• 


Thus, it is readily seen that 



o 


(4-27) 


(4-28) 


Since the compound nucleus can decay back into the incident 
channel, the summation in Equation 4-28 has to be over the incident and 
all the outgoing channels. Thus, combining Equations 4-24 and 4-28 yields 


o 


a|i 



(4-29) 


Equation 4-29 is essentially the Hauser-Fesbach formula for the case of 
spinleas particles. The generalization of Equation ‘,-29 requires the 
Introduction of the spins of the channels involved. Let i denote the 
spin of the Incident particle, 1 the spin of the target nucleus and 
j = I + 1 the incident channel spin, and j ' s I' + i 1 the corresponding 
quantity in relation to the exit channel 8, where i' denotes the spin 
of the emitted particle and 1' that of the residual nucleus. Let i and 
Z 1 denote the orbital angular momentum in the incident and outgoing 
channel, respectively. 

The entrance channel a is specified by i, I, j, and 1, which 
combine to give a resultant: singular momentum J. With these additional 
specifications on spins, the transmission coefficient T a is written as 
T a £j . The probability for the formation of a compound nucleus of spin 
J from the channel o is then given by 

a (CN)J = „ *2 (21 + 1) |C J J £ | 2 » (4-30) 

a a al j 1 mmo 1 


where 

mmo 


are the well-known Clebsch-Cordan coefficients. 
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The details of this equation with justifications of the various 
terms entering in the calculations are worked out by several groups and in 
particular by Blatt and Beidenhorn (Ref. 12). The most Important terms 
that enter into Equation 4-30 are T^^ , which are the transmission 
coefficients. These transmission coefficients are calculated by the choice 
of an appropriate optical potential. Once the compound nucleus is formed 
^nd the related cross section is established through Equation 4-30, it 
only remains to write down the cross section for the decay of the 
compound nucleus through the exit channel. Hauser and Fesbach (Ref. 13) 
and Hodgson (Ref. 14) have discussed the details of the cross section and 
have given 





U ± 1 

(21 + 1 ) (21 + 1 ) 


Aj (ji, j * , f. 1 ; 


6 ) 


T 

q*j gr,r 



» 


(4-31) 


where the Aj coefficients are products of the Clebsch-Gordon and Racah 
coefficients and are tabulated by Blatt and Beidenhorn. 

Without going intc the details of the theory, which is rather well 
formulated, it is of interest to mention the role of optical potentials 
that enter into the calculation of the transmission coefficient T^ . In 
any calculation using this formalism, the accuracy of the results would 
depend a great deal on the appropriateness of the optical well parameters. 

4.2.4 Optical Potentials for Charged-Particle Reactions 

Ever since Serber (Ref. 15) first proposed the optical model 
for the nuclear potential as seen by an approaching nucleon, this model 
has received more attention than any other form. The resemblance to a 
complex refractive index that appropriately explains reflection and 
refraction can very well explain the physical process of elastic collisions, 
which are similar to reflection, and to compound nonelastic collisions, 
which are similar to the phenomenon of refraction. There are several 
forms of optical potentials. In the analysis in this study, the following 
general form is used: 
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V(r) ■ U((r) + lWg(r) + (Uh + IW») 


MiM l • a . 

dr 


It is also found that the actual form of the potential, as well 
as the variable parameters occurring in Equation 4-31, are dependent 
on the charge and energy of the projectile. For the case of proton 
reactions with light nuclei and for protons of medium energy from 9 to 
22 MeV, Perry (Ref. 16) has given the following values for the various 
parameters occurring in Equation 4-32: 

U - 53*3 - 0-55E - 27 5-— 3. + 0-4 ^ 


0, 5 for E > 17 


U 8 - 7-5 for E < 17. 

In Equation 4-32, the tern corresponding to U 8 comes out to be 
„ , . / 11 \ 2 „ 1 df (r) r - 

v - (r) ■ (sHc) "s 7 d— r • 

This is called the spin orbit term, where / -™V “ 2 • fm 2 is used in the 

\ m |[f' / 

calculations, L and o have the usual meaning respectively, and f(r) is a 
form factor. 


Some typical values of the various energy-dependent parameters 
used on the case of proton reactions are given in Table 4-1. 

The mechanism and the components of the various possible types of 
channels in the compound-nucleus formalism are shown in Figures 4-1, 4-2, 
and 4-3. The cross sections for various types of reactions referenced 
in this discussion as well as in most literature pertaining to these 
reactiot.a are shown in Figure 4-4. An ideal formalism must be able to 
calculate each of these components as a function of the projectile 
energy and spin and the channel parameters. However, the mechanism 
to be used for the direct reaction o(di) is definitely different from. 
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TABLE 4-1. PARAMETERS FOUND IN THE ANALYSIS OF PROTON DATA WITH DIFFERENT POTENTIAL FORMS 
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for example, the Hauser-Fesbach formalism, which will be used in this 
discussion for the evaluation of inelastic cross sections. 


4.2.5 Hauser-Fesbach Calculation* 


Based on the formalism developed in the previous section and 
adopting ABACUS 11 (Ref. 17), the inelastic cross sections are calculated 
for 8®^* 7^^* anc * 6 C ^ • However, the following modification has been 
made in the existing ABACUS II. The coulomb subroutine has been replaced 
by the more recent Manchester Code (Ref. 18). The Coulomb wave functions 
Fi(n,p) and Gj^(n,p) and their respective first derivations are determined 
from the solution to 


di_ 

d p 


U* 


In _ iq + i) 

p " p 2 


o , 


where 


P 



r 2 \V 2 

& ) 


n 


Z -l.. Z 2-- e 2 

1i ? k 


where Zje is the charge of the projectile and Z 2 e that of the target, 
k is the wave number, and E is the cent r-of-mass energy. The program 
for the coulomb wave function can be incorporated into the ABACUS II 
or may be run separately and the relevant data put into die calculation 
to be made for charged particle leactions. 

The intrinsic parameters of the target nucleus for oxygen 16 
are shown in Figure 4-5. The various parameters that enter into the 
optical potential well used in the HF calculation are given in Table 4-1. 
The energy dependence of these parameters is clear from Table 4-1. 

Table 4-2 gives a comparison of these parameters at a fixed proton 
energy of 16*3 MeV. ABACUS II in conjunction with the Manchester 
Code strictly calculates only the (p,p') reaction cross sections. There 
is no provision in the formalism developed to handle any nonelastic 
channels leading to charge exchange or any other type of inelastic process. 
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TABLE 4-2. PROTON INELASTIC CROSS SECTION FOR g 0 16 


\ Ep 

0 (p.p') in mb 


11 

12 

13 

14 

15 

6-06 

16-7 

unr 

29-8 

47 3 

56 

6 14 

12.3 

16-4 

25-6 

31-6 

36-5 

6-92 

12-2 

11-2 

17-3 

22-4 

25-8 

7*12 

9-7 

9-9 

12-6 

16-2 

18-8 

8-88 

6-2 

6-8 

10-9 

16-2 

16-5 

9-58 

3-3 

3-9 

6-8 

9-7 

10-3 

9-84 

9-6 

9-9 

16-6 

21-4 

19-7 

10-36 

12-2 

13-1 

19.9 

27-3 

30-1 

10-94 

0-0 

0-2 

1-0 

1-5 

2-0 

11-07 

6-8 

6-9 

9-6 

16-2 

12-9 


The cross sections to excite each of the states are given in Table 4-3 
as a function of Ep, where Ep denotes the proton energy. 

i a 12 

Tables 4-4 and 4-5 give similar calculations for ^N 10 and gC ; 
Figures 4-6 and 4-7 show the respective transit ionr. of interest. 

4.2.6 Discussion o £ R esults 

The results of the calculations carried out by adopting the 

modified ABACUS Code with the introduction of the Coulomb routine through 

1 7 

the Manchester Code are given in Tables 4-2, 4-3, and 4-4 for ^C , 

^N^’, and gO^’ nuclei and for proton energies from 11 to 15 MeV. The 
parameters of the excited states entering the computation are shown in 
the respective local schemes of these nuclei in Figures 4-5, 4-6, and 4-7, 
respectively. The calculation is clearly dependent on the optical 
potential parameters chosen. Bechet ti and Greenlcss (kef. 21) have 
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TABLE 4-3. PROTON INELASTIC CROSS SECTION FOR 7 N™ 


M 



< 

j(P.P') 

In mb 



11 

12 

13 

14 

15 

16-5* 

16-6** 

2-31 

8*8 

6*9 

5*7 

7*9 

6*6 

2*9 

3*2 

3-94 

18-2 

23*4 

18*9 

16*9 

22*3 

17*5 

18*5 

4-91 

5*3 

4*3 

6*1 

4*6 

4*4 

3*3 

1-8 

5*1 

38*26 

29*0 

32*1 

29*6 

26*9 

16*0 

19*0 

5-69 

18*2 

16*6 

11*96 

16*9 

16*0 

8*8 

8*7 

5-83 

56*9 

43*1 

32*7 

45*38 

40*1 

19*6 

19*2 

6-23 

10*28 

13*7 

11*59 

16*7 

13*8 

5*8 

1*7 

6-60 

10*1 

16*8 

11*9 

16*9 

15*3 

6*0 

5*9 

7-03 

19*4 

15*4 

23*2 

17*2 

26*9 

23*6 

27*9 

7-4 

10*2 

8-4 

11.7 

9*4 

9*6 



7-6 

37*7 

28*7 

32*7 

29*4 

27*4 



7-96 

10*65 

16*1 

13*4 

17*4 

16*4 



8-06 

5*7 

8*4 

7*2 

9*3 

8*9 



8-62 

1 *8 

2*6 

2*3 

3*0 

2*8 




♦From Reference 19 
♦♦From Reference 20 
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TABLE 4-4. PROTON INELASTIC CROSS SECTION FOR 6 C 12 


\ Ep 
EX(MeV) 
(MeV)\ 

o(p.p') In mb 

11 

12 

13 

14 

15 

4-433 

31-8 

51-2 

59-6 

64-7 


7-65 

6-6 

7-0 

8-2 

10-1 


9-63 

20-0 

25-0 

27-6 

31-7 


10-10 

1-6 

2-2 

4-2 

4-6 


12-73 

3-9 

4-9 

5-9 

6-5 



TABLE 4-5. OPTICAL MODEL PARAMETERS USED IN THE CURRENT CALCULATION 


OPTICAL MODEL 
PARAMETERS* 

mm 

m 14 

7 n 

J6 

8° 

VRC 


50-4 

51-62 

VIM 

19-9 

3-0 

3-0 

R1 

2-747 

3-05 

1-25 

A1 

0-4 

0-68 

0-58 

RC 

2-74 

0- 

1 -25 

R2 

2-74 

3-05 

1-25 

A2 

0-25 

0-7 

0-6 

SO-REAL 

6-0 

7-0 

6-9 


*The various symbols used correspond to the 
description in ABACUS II 


« 
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carried out extensive analyses of potential parameters for proton energies 
less than 50 MeV. The most complete analyses h:.ve been carried out b> 

Perry (Ref. 16) for proton elastic scattering from 9 to 22 MeV. A summary 
of these analyses Is shown in Table 4-1, which is adopted from Hodgson 
(Ref. 14). The values of the various optical model parameters defined 
in ABACUS are given in Table 4-2. 

The calculated values of the Inelastic cross sections lor each 
of the states of 7 N*^ at a proton energy of 14 MeV, as given In Table 4-3, 
compare favorably with reported results. lowever, there Is some degree 
of disagreement because the lnelustlc processes that open up due to other 
types of reactions are not considered In the formalism. Therefore, 
at best this technique can calculate only the inelastic cross sections 
for the case in which such particle emission cross sections are negligible. 

in the case of light nuclei of Interest here, for energies of about 
14 MeV, other inelastic channels can become available. Therefore, It is 
found to be necessary to incorporate a theory for handling such situations. 
Such a forma 1 Ism is discussed in the next section. 

4.2.7 Coupled -Channe 1 Formal Ism 

The application of compound-nucleus formalism witli no regard for 
exit chunnels leading to emission of different particles Is, at best, 
rather a crude approximation when the projectile energy is of the order 
of or higher than the observed Q value for a particular reaction. However, 
as mentioned earlier, if such reactions are to be taken into considera- 
tion, other theoretical formalism should be developed. 'Hie vaii.it ion 
of cross section with the energy of the projectile as well as reported 
shape of the angular distributions lead to the conclusion that, at 
energies of interest here, the contribution from the compound state 
to other channels is not negligible. Thus, for example, in Table 4-3 
the calculated value and the measured value of inelastic cross sections 
of 7 N 1/ * at an incident proton energy of 16.5 MeV are not in agreement. 

The calculation employing the Hauaer-Fesbach formalism will give the 
compound inelastic cross section only. The direct inelastic cross section 
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leading to Inelastic channels, such ns (p,ut) reactions, will have to be 
handled by adopting the macroscopic and microscopic coupled-channel cal- 
culation as developed in the Oregon State ” < -..orsity coupled-channel code 
(Kef. 22). 

4.2.8 Coupled Channel Theory 

In a simple situation one can assume that a spinless particle 
reacts with a nucleus and raises the nucleus to a series of spinless excited 
states . 


Let the kinetic energy operator corresponding to the incident 
particle be presented by 


T 


Hi 

2m 


V 2 

r 


» 


the nuclear Hamiltonian be given by II (fj , and the interaction potential 
between the projectile and the nucleus be represented by V(r, (,) . The 
Schrodinger equation pertaining to this situation can be written as 

[T - V(r, 5 ) + Ha)] ¥(r, 6) - E'Kr, O . (4-33) 


Tile nuclear states must be eigen states of the nuclear Hamiltonian so 
that these states, XuCO » should satisfy the relation 

HU) x a (0 * e a Xa(^) • (4-34) 

The eigen functions, XaCOt must form a complete orthonormal set, and hence 
the wavefunctions can be expanded in terms of this orthonormal set and 
can be expressed as 

nr, O « E T (r) 7(0 (4-35) 

a a ° 

It now becomes possible to make use of Equations 4-34 and 4—35 
in Equation 4-33. By substituting Equations 4-34 and 4-35 in Equation 4-33 
and after multiplying on the left by x *(C), the resulting function is 

integrated over the entire nuclear coordinates to obtain 
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**(*;) [T - V(r, O + H<0]£ ^(r) X «U> d£ 
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V a (r) x a U) d£ 


(4-36) 
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H'a(r) X a (5) 
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= / *2 V I*. 
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(4-37) 


In Equation 4-37 
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aa 


(r) 



X '*(?,) V(r , O x’(0 df, 

£1 d 
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The usual technique of expanding into partial waves can be used 
at this stage, and thus one writes 

V (r) = S’ Y* (0,*) . (4-38) 

Thi6 expansion is substituted into Equation 4-37, then multiplied on the 

left by y!' 1 *, and then integrated over 0 and $ to obtain 
L 


y*“*(e.*) (T - e 0 ) £ y“ (e,*) d ft 


N 


/ M* r' U,(r) m 

*" <“•« 2 ~ f h < 0 '»> da • 


a’=l 

LM 


(4-39) 


Here 


E a = E * e a 


*fik 2 

2m 


Let 
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(4-40) 


Then Equation 4-39 can be written as 


iik 2 

T " 2m” ” V 


U a (r) = w aa -( r ) U a (r) 


This equation can be written in the form 


ii 2 d 2 Ti 2 ?. (? + 1) „ „ .. . r F 

2^ dF ' 2^P “ V 0pr V S.O V coul + L E n 


U a (r) 


■ I V i„t U a (r > ' 


(4-41) 
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Equation 4-41 couples all the elastic channels with the inelastic channels. 

The equation can be numerically solved adapting the coupled- 
channel code developed by Oregon State University. It is of interest to 
pursue this aspect of the problem to obtain a meaningful result fo** 
the various inelastic cross sections, particularly those leading to the 
particle emission channels. However, it is expected that such a study 
would consume several months of effort. 

4.3 CONCLUSIONS AND RECOMMENDED PLAN FOR FURTHER STUDY 

The coupled-channel code as developed by Oregon State University 
is a rather extensive analytical technique. The original code was 
developed for the CDC 6600 assembly. It has been possible in this 
effort to study the applicability of this formalism as an improvement 
on the possible inelastic cross section calculations. It is strongly 
felt that ABACUS II in conjunction with the Manchester Coulomb Code and 
the coupled-channel code can lead to a more extensive and complete 
analysis of the charged-particle reaction in the case of light nuclei. 

The coupled-channel code has been adapted for the Univac 1108 assembly 
and is ready to be put into operation should such an investigation be 
considered of interest. 
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